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Abstract 

The entropy per baryon is studied in the strong couphng large di- 
mension d hmit of lattice QCD with staggered fermions. The partition 
function is calculated for non-zero chemical potential and temperature 
using the 1/d expansion. It is found that the entropy per baryon ratio 
is almost continuous across the transition from the the quark-gluon 
to the hadronic phase. The relevance of these results for heavy ion 
collisions is discussed. 

Relativistic heavy ion collisions test the thermodynamic properties of 
QCD. It is well known that lattice gauge theory is the most suitable method 
for studying these properties theoretically, however Monte Carlo simulations 
at non-zero chemical potential /i are notoriously difficult ||, 0], and it is 
necessary to investigate separately analytical methods based on the strong- 
coupling, 1/d expansions of QCD at both T ^ and /i 7^ 0. This is the 
subject matter of the present investigation. This method has been success- 
fully applied to the calculation of the hadronic spectrum [|| and leads to a 
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quantitative description of the phase diagram in the temperature and chem- 
ical potential plane |l|, |], H, ^ |^. It yields a second order phase transition 
at vanishing fi and predicts a first order transition for any non-zero value 
of fi. However, unlike the weak coupling regime, where the dimensionless 
lattice parameters are related to the lattice scale via the renormalization 
group equation, such a relation does not exist in the strong coupling regime. 
Consequently, a physical interpretation of the strong coupling results is not 
straightforward. In this paper we propose a natural solution to this problem. 
We calculate the thermodynamic properties of QCD and contrast the results 
to those obtained in the hadronic gas model. The difference in results is 
very remarkable. In the ideal gas approximation the entropy to baryon ratio 
(S/B) is almost always larger in a quark-gluon plasma than in a hadronic 
resonance gas, if we compare them at the same temperature and chemical 
potentially, |rT|. Our results provide no support for the contention that 
experimental results on the charge asymmetry ratio show evidence for quark- 



gluon plasma formation |12 



As is usual in these considerations, we work on an asymmetric lattice with 
Nt points in the time and Ng in the d space directions with an anisotropy 7. 
The finite-temperature behavior can be studied in the infinite volume limit 
{Ng — > 00) as a function of 7 and Nt. In the naive continuum limit 7 can 
be interpreted as a ratio of the space and time lattice spacing, as/ at. As 
we shall shortly demonstrate, this naive relation does not hold in the strong 
coupling regime. We shall propose instead a new relation based on general 
physical considerations. 

We start from the effective partition function derived at infinite coupling 
and in the large-c? limit 



Z^= [ [dX] exp{--^ ^(A,. - m)V-y\Xy - m)} [] / dU,detD{Us), (1) 

x,y X 

with 

1 

^^,y = ^ ^ i^y,x+k + ^y,x-k)- (2) 



k=l 
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D{U) is an NtNc x A^^A^^ matrix 
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with Ui = 2(Aj + m)/7. In the final expressions the parameters and at 
each space-time point x are set equal to the constants a = e'^ and b = e~'^ 
which are related to the chemical potential. The explicit x-dependence of 
the Oj, bi and Aj in has been omitted. The integration over dU can be 
performed by exploiting the appropriate SU(A'"c) group integral. Details are 
given in appendices of refs. P, |T3 |. 



The 1/d expansion begins with a saddle-point approximation of the A- 
integral, or equivalently by minimizing the mean-field free energy 



d 



(4) 



with Zq = J dUdetD{U) being a "zero- dimensional" partition function. We 
improve infinite-coupling mean-field calculations by replacing the free energy 
(i by 



1=1 



(5) 



with corrections SFi of order or 1/d 10, |^ . 

Provided we know how to fix the scale and at we can relate the physical 
thermodynamical observables to the corresponding lattice quantities derived 
from the free energy (|,^ 
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where 



F ^ I OF 

OF 



latt nr Pvac i ^latt t\t t\t ' 



^latt OAT ^vac 5 (7) 

Slatt = Nt{Platt + ^latt- l^Njliatt)- 

Pyac and Syac are the vacuum pressure and energy density calculated in the 
limit Nt oo. ub and s refer to the baryon density and the entropy density 
respectively. The entropy per baryon ratio S/B is then given by s/hb- 

At /i = the mean-field calculation predicts a second-order chiral phase 
transition. Including 1/d and l/g"^ corrections yields the critical anisotropy 

where the critical anisotropy at infinite coupling takes on the value 

_ d{K + l){N, + 2) 
^0- 6(iV, + 3) 

The anisotropy 7 is a measure of the temperature. However, unlike the 
weak coupling regime, where 7 can directly be related to the spatial and 
temporal lattice spacings |T^ allowing one to establish the connection with 
the temperature T = 'j/Nfas, such a direct relation does not exist in the 
strong coupling limit. In other words, the physical anisotropy of the lattice 
as/ at is not necessarily equal to 7. If one insists on using this relation in 
the strong coupling regime the critical temperature, owing to (||,|D, would 
strongly depend on Nt which is physically unacceptable. We must, therefore, 
alter the naive relation between 7 and the physical anisotropy (or tempera- 
ture), in such a way as to maintain the critical temperature independent of 
unphysical lattice parameters, such as A^^^. Quite generally we may set 

(10) 

where the function / can be determined from the requirements: 
a) Tc should not depend on Nt, 
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b) at 7 = 1 we have at{l,g) = as- 
From the first requirement it follows that the Nt dependence of /(7c, g) must 
be of the form Nt(p{g). Furthermore, in the limit g ^ oo owing to (|^) we 
must have /(7o,oo) = 7q up to a constant which may be absorbed in the 
lattice scale a^. Thus /(7c, g) = 'yQf{g), with ip being an arbitrary function of 
g. From the functional relationship between 7c and 70 in (H) it immediately 
follows 

/h.j)=7= + ^(2d-^). (11) 

up to an arbitrary multiplicative function of g. The second requirement then 
implies 

ath,9) = as^r^y (12) 

We note that the correction to the critical anisotropy does not depend 
on Nf although 7q itself is proportional to Nf. The temperature can now be 
defined in the usual way T = 1 / Ntat{'y , g) . Thus, for example, at infinite 
coupling we have T = 'y^/N^as in contrast to the standard T = 'y/Nfag. The 
physical chemical potential is related to the dimensionless lattice parameter 
H as usual Uphys = fJ-atij, g). 

The scale in (|T2|) remains completely arbitrary. Unlike in the weak 
coupling regime, where the scale can be related to the coupling constant via 
the renormalization group equation, such a direct relation does not exist in 
the strong coupling regime. One can, nevertheless, fix the scale by fixing the 
value of some physical observable, such as or /ic, at a chosen value of g. 

In Fig |1| we show the critical temperature as a function of 6/g^ for a fixed 
l^phys- The reason for a slight Nt dependence is that the numerical evaluation 
of the transition temperature includes higher powers of l/g"^ which are not 
present in the analytic formula (§) used to derive (|Tl]). The scale is fixed by 
requiring that the strong coupling Tc at /i = and 6/g^ = 5 coincides with 



the Monte Carlo value T,. = 130 MeV 15 



The critical entropy per baryon in the hadronic phase increases rapidly 
with decreasing g (Fig If we fix the coupling at the point where the 
strong-coupling curve crosses the experimental region (hatched area) |]16| we 



can plot the entropy per baryon as a function of temperature on both sides 
of the transition point and compare this with the quark-gluon plasma and 
hadronic gas model (Fig H). In contrast to this simple phenomenological 
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model, strong coupling lattice thermodynamics predicts a continuous s/ub 
across the phase transition. In the infinite coupling limit (6/(7^ = 0) the 
entropy per baryon in the chiral-broken (hadronic) phase is greater than in 
the symmetric (quark-gluon) phase. 



Preliminary data from the emulsion experiment EMU05 cited in [|T6l show 
that in S-Pb collisions at 200 GeV the charge asymmetry ratio 

(13) 

is approximately 0.085±0.01 in the central rapidity region. To relate this 
to the entropy per baryon ratio one can proceed as follows. Assuming that 
this region is symmetric in isospin (for S-S collisions this would be exact), 
the numerator in (0) equals the baryon number divided by 2, assuming 
furthermore that the denominator is approximately 2/3 of the total particle 
number (neutral particles are not included), the ratio becomes : 

« ^ (14) 

To relate this to the entropy we use as a guess that the entropy per particle 
is 4 (as is the case of a massless Boltzmann gas). We thus conclude that 

S/B ^ 36 ± 5 (15) 

This estimate is compatible with the more detailed analysis presented in[p!6|. 



T7| . This result was used above in Fig. 2. 

Contrary to the predictions of quark-gluon plasma - hadronic gas models 
we find in the strong coupling regime that the entropy per baryon is almost 
continuous across the phase transition. Moreover, including the strong cou- 
pling corrections increases the entropy per baryon content in the hadronic 
phase near the phase transition to the order close to the estimates obtained 
from the measured charge asymmetry ratio. 

In conclusion, we have calculated the thermodynamic properties of QCD 
in the strong coupling large d limit at finite temperature and chemical po- 
tential. We have introduced a new relation between the lattice anisotropy 
parameter 7 and the physical anisotropy a^/aj allowing us to establish a 
physical interpretation of the dimensionless parameters 7 and /x in terms of 
temperature and the chemical potential. Our results show that the entropy 
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per baryon ratio S/ B is almost continuous across the phase transition point. 
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Figure 1: Critical temperature at fixed jj, — and — 250 MeV for Nt — 4 
and 6. 



Figure 2: Critical entropy per baryon in hadronic phase at fixed /j, — 250 
MeV for iV^ = 4 and 6. 



Figure 3: Entropy per baryon calculated in strong coupling QCD compared 
to quark-gluon plasma (QGP) and hadronic gas model at fixed /x=250 MeV. 
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